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Abstract 

Electroweak baryogenesis provides a very attractive scenario to explain the ori- 
gin of the baryon asymmetry. The mechanism of electroweak baryogenesis makes 
use of the baryon number anomaly and relies on physics that can be tested experi- 
mentally. It is today understood that, if the Higgs mass is not larger than 120 GeV, 
this mechanism may be effective within supersymmetric extensions of the Standard 
Model. In this work, we reconsider the question of baryon number generation at 
the electroweak phase transition within the context of the minimal supersymmetric 
extension of the Standard Model. We derive the relevant diffusion equations, give a 
consistent definition of the sources, and compare our results with those appearing 
in the recent literature on this subject. 
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1 Introduction 



Electroweak Baryogenesis PP provides a predictive framework for the computation of 
the baryon asymmetry of the Universe [2]. Perhaps the most attractive feature of this 
mechanism is that it rehes on anomalous baryon number violation processes which are 
present in the Standard Model At temperatures far above the electroweak phase 
transition scale, these anomalous processes are unsuppressed and, in the absence of any 
B — L asymmetry, they lead to the erasure of any baryon or lepton number generated at 
high energy scales j3] . These baryon number violation processes are, instead, exponentially 
suppressed in the electroweak symmetry broken phase, at temperatures far below the 
electroweak phase transition [H]. At the electroweak phase transition, non-equilibrium 
processes may generate a non-vanishing baryon number which may efficiently diffuse into 
the broken phase jH] • The generated baryon number depends on the CP- violating currents 
present in the model. The mechanism of electroweak baryogenesis may become effective if 
the CP-violating sources are strong and, at the electroweak phase transition temperature, 
the baryon number violation processes in the broken phase are sufficiently suppressed, 
leading to a baryon number density in the broken phase consistent with observations. 
This, in turn, demands a strongly first order electroweak phase transition. 

It has been long ago realized that in the Standard Model the CP violating sources are 
too weak to lead to an acceptable baryon number density [3 E]. Moreover, even if the 
sources were strong enough to lead to a reasonable baryon number generation, the phase 
transition is weakly first order, leading to a strong exponential suppression of the baryon 
number generated in the broken phase jH]. A minimal supersymmetric extension of the 
Standard Model (MSSM), instead, has all the necessary ingredients to improve on both 
problems First, there are additional sources of CP- violation, associated with the 

CP-violating phases of the supersymmetry breaking parameters. Second, in the presence 
of a light stop, the phase transition may become strong enough to allow the preservation 
of the baryon number generated at the electroweak symmetry breaking scale [TT]-|28j. 

Although there is general agreement on both the existence and nature of the new 
sources and also on the parameter space leading to a strongly first order phase transition, 
the results regarding the strength and specific form of the CP-violating sources are still 
controversial. Most groups working on the subject have found that, for values of the 
CP-odd Higgs mass of the order of the weak scale, the sources are dominated by a term 
proportional to eijHid'^Hj, where Hi denote the expectation value of the neutral compo- 
nents of the two Higgs doublets present at low energies [29^-^39]. Recently, however, it has 
been claimed that this contribution is absent, leading to a suppressed result for the baryon 
asymmetry contribution within the MSSM proportional to Hid^H2 + (1^2) pT| I42j . 

In this article, we proceed to perform a critical re-analysis of the diffusion equations 
leading to the baryon asymmetry and of the sources for those diffusion equations. In 
particular, we demonstrate that the sources can be expressed in terms of appropriate 
differential operators times the CP-violating currents we computed in a very recent work. 
We find no suppression of the dominant sources proportional to eijHid^Hj. Indeed, we 
show that the suppression claimed the authors of Refs. [HI US] is an artifact of the 
approximation they used in order to compute the currents. 

The organization of the article is as follows. In section 2 we provide the derivation of 
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the diffusion equations. In section 3 we review the derivation of the CP-violating currents 
and we show exphcitly where our results differ from the ones of Ref. [H]. In section 4 we 
discuss the results for the baryon asymmetry within the MSSM and section 5 contains our 
conclusions. Finally some useful formulae concerning the chargino sector of the MSSM 
are summarized in appendix A. 



2 Diffusion equations 

We will start considering a system of particles propagating in a non-trivial CP-violating 
background localized in the bubble wall, where the bubble is expanding with a speed 
in the plasma frame. The presence of the CP-violating background and the particle 
number changing reactions perturb the distribution functions for particles and antipar- 
ticles around the equilibrium. The corresponding distribution functions fi satisfy the 
Boltzmann equations, 

v^'d,U + FtV,U = C,[f] (2.1) 

where we are using the notation 9^ = d/dx^, = d/dp^, v'^ = dx^/dt = p^/E is the 
four- velocity and = dp^ /dt is the force generated by the non-trivial background. The 
term on the right-hand side Cj = {dfi/dt)^^^^ encodes the effect on the distribution of the 
particle number changing reactions and elastic collisions. 

In the bubble-wall frame and neglecting the curvature of the bubble wall the distribu- 
tion functions can be written as, 

rP _ 9{Pz, z) \ SfP (9 9\ 

where /if = /if {z) is the chemical potential, E and pz are the energy and momentum of 
the particle, the momentum-dependent part df^ = 6ff{E,pz;z) describes departure of 
the system from kinetic equilibrium and we have introduced a slowly varying CP-even 
function of Pz and z, g{pz,z), which tends to 1 in the absence of forces. Considering the 
effective distribution functions fi = ff — ff, describing the difference between particles 
and anti-particles of a given species, and applying ()2.2|1 to Eq. ()2.1|1 gives 



El 
E 



{F^ + v.F:)-g{pz,z)^ = a (2.3) 



where Fj^ is the CP-violating component of the force, we have linearized with respect to 
the small perturbations, /ij, 6fi, in the spirit of the Chapman- Enskog first order approx- 
imation /o is the equilibrium distribution function in the plasma frame, and the (') 
denotes derivative with respect to z. In this process, due to the assumed smoothness of 
the CP-even function g{pz,z), we have ignored terms proportional to derivatives of this 
function. 

Following Ref. jH], we rewrite the above equation in the plasma frame by performing 
a Galilean transformation on the velocity Vz to Vz — foj, 



Pz 

E 



-fiig{pz,z) -Q^ + ^fi 



+ {F^ + v.F:)g{pz,z)^ = a (2.4) 
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where g{pz, z) is assumed to be an even function of momenta in this frame. 
We now define 





(2.5) 



(2vr) 



and using the decomposition 



f^i^g{Pz,z) + 6fi 



(2.6) 



and the definition of the current 



(2.7) 



we can relate the CP-violating current with the chemical potential and the function 6fi 

as, 



where, again, we have neglected terms proportional to derivatives of g{pz, z). The function 
6fi is such that (Sfi) = and the number density is entirely provided by the chemical 
potential /Xj. In this way we can choose, with all generality, the function 6fi as odd with 
respect to Vz and therefore satisfying the property {v1^5fi) = 0. 

We now multiply Eq. ()2.4|1 by v'^ = {l,pz/E) and integrating over momenta we get 
the two equations. 



where Tj is the total interaction rate and the diffusion coefficient Di is defined as 



and we have again ignored the smooth momentum dependence of g{pz, z) in the numerator 
and denominator integral functions. Observe that the ratio ()2.1H1 remains constant up 
to derivatives of g'{pz, z) which we have consistently ignored in our treatment. 

The terms on the right-hand sides of Eqs. ()2.9|) and ()2.10|) can be decomposed, using 
the decomposition ()2.6|) as 




(2.8) 





(2.11) 



(2.12) 
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where Cj and Cf are linear in /ij and 6fi, respectively. In fact (C/[/i]) is entirely provided 
by inelastic collisions 

71 ■ 

(C,^H>^-r,,-l (2.13) 

where kj = 2 {kj = 1) for bosons (fermions) are statistical factors, Tij (j = ji,j2,...) 
is the averaged interaction rate for the inelastic reaction channel and rij/kj implies 
a sign over all particles participating in the reaction with appropriate sign assignments: 
For an inelastic reaction corresponding to the decay i jij2 with decay width Tij^j^, its 
contribution to the right-hand side of ()2.1Hj) would be —Tij^j^{ni/ki — rij^/kj^ — rij^/kj^). 
Furthermore {Cf[6f]) = from the oddness of 5/,. On the other hand, {vzCHfj]) = 0. 

Within the present framework, the CP-violating forces are provided by the interaction 
of the different fields with the Higgs profiles, which vary along the bubble wall. These CP- 
violating forces induce CP- violating currents for the different fields, which lead to {rti^jf) 
after elastic and inelastic interactions with the other fields present in the plasma. The 
resulting currents have then two components. The first one, {n^^\ the background 
component, is the one that would be obtained in the presence of only the interaction 
with the background Higgs field, and the second one, the "collision" component comes 
from the interaction with the plasma fields and lead to the diffusion process. The integral 
{vz C^[5f]) is contributed by all collisions and dominated by the elastic ones. Correspond- 
ingly with the two components, background and collision, of the currents there can be 
defined two 5/ components. 

The part of 5f that is induced by the interactions with the plasma fields, the collision 
component, is governed by the same interactions as the ones governing Ci and is therefore 
relevant in the same regime of momenta in which the collision term becomes important. 
On the other hand the 5f component coming from the interaction with the background 
fields, 5f^^\ is not correlated with the collision terms and leads naturally to a negligible 
contribution to Cf[6f]y Therefore, considering that the rate of interactions is ap- 
proximately constant in the momentum regime in which the collision contribution to the 
currents become relevant, and negligible anywhere else, we get 

(|cf[5/])^rf (jf-jf)^). (2.14) 

Using now ()2.13|) and ()2.14|) in ()2.9p and ()2.10|) we obtain the diffusion equations, 

-vX + {jD' + (^{F^ + v..F:)g{pz,z)^^=-T,,f^ (2.15) 

D^K - UD' + (I + -^^n 9ip^^ -)||) =rf (jT - J^) ■ (2.16) 

A current j^^'' ^ = (^n\^\ ^ j in the presence of a CP- violating Higgs background (or 
CP-violating force F^) but where no interactions with the plasma fields were considered 
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was computed in Ref . [S^ . It should satisfy the set of Boltzmann equations 



-v^ {nf')'+{r'')'+(^{F^ + v^F:)gip,,z)^'j =0 (2.17) 
Arf (nf^)'-.. (jf ^0'+ (I i^' + -^F:)gip.,z)^'^ =0 (2.18) 
We can see from these equations, and in particular from ()2.17|) that the roles played by 
()2.17p from ()2.15p one obtains an equation 



the CP- violating forces and currents {nl , ) are equivalent. In fact by subtracting 



-v^ (n, - nS^))' + (,T - ,f )^)' + r,,^ = (2.19) 

where the force terms have been replaced by the divergence of the CP- violating current. 
By subtracting now ()2.18|) from ()2.16|) we obtain, 

- 3^') = D. [(n, - „f >)' - ^ („, - „f' )" + 0(vl)] ^ A - ,4"')' (2.20) 

where the last approximation comes from the requirement Vi^/TjL^ <C 1 (where is 
the wall thickness), which is necessary for the validity of the derivative expansion we will 
use in the calculation of the CP-violating currents in section El Eq. ()2.20j) represents the 
well-known empirical law of diffusion or Pick's law. 

Replacing now Eq. ()2.20j) into ()2.19|) we obtain the final expression for the diffusion 
equation as, 

-v^n[ + D.n'l + P,,^ = (2.21) 

Kj 

where the source is given by 

S. = A [nT'X - V. inTA' ^ A fnf ))" (2.22) 



and the last approximation comes from the fact that Di ^ L^v^)- Notice that Eq. 1)2.211) 
has the correct boundary conditions since in the absence of inelastic reactions (for Pj^ = 0) 
it provides the trivial solution = n] as required. 

We can compare this result with the source obtained in Ref. jHj. These authors work 
in a WKB approximation with quasi-particles. Within their semiclassical approximation 
the energy is a (constant) label of WKB-states and therefore the time component of the 
force = dE/dt vanishes, while they have assumed that g{p;^,z) = 1. Moreover their 
source, Sf'^^, is given by IH], 

^P" = -^^(|-^^^||y (2-23) 

which appears in Eq. ()2.18)1 above and can therefore be related to n^-^"* and ^ as they 
also appear in ()2.18p . Although we have obtained the sources in a way independent 
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of the form of the CP- violating forces and for vanishing we obtain a relation 
between the sources Si and Sf"^^. Indeed, using Eq. ()2.18j) . we obtain 

Sf«==A(«f')"-^(.f")'-S<- 

where the last approximation relies on the relation DiVj ^ 1. The latter relation holds 
for the case of charginos in the MSSM, analyzed in section IHl 



3 The sources 

In this section we will make contact between the source in ()2.22|) . and in particular the 
current j^^^^, and the Green function and will apply the formalism to the sector of 
charginos in the MSSM. 

First consider a chiral fermion (say a right-handed one) in the presence of the non- 
trivial background. Its Green function S^^\x,y) can be considered as function of the 
center-of-mass coordinate z = {x + y)/2 and the relative coordinate r = x — y, 

S^^\z + r/2,z-r/2). 

Since we will make a gradient expansion with respect to the coordinate z we perform a 
Fourier transformation with respect to the relative coordinate r as 

SiB){z;p)= /rfV e'P" 5(^)(z + r/2,2-r/2) (3.1) 



which is called the Wigner representation of the Green function S^^\ 
By making the general decomposition 

= (T^5(^)(^;p) (3-2) 

one can define the corresponding background current j ^ as 

and making the inverse Wigner transformation (jH.lj) it can be cast as 
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-limTrcT^5(^)(z + r/2,2-r/2) (3.4) 



2 r^O 



The current can be written in an alternative equivalent form, that makes connection 
with the discussion in the previous section. 



] 



Under small perturbations the relevant component of the integrand is approximately 
dominated by 5 functions imposing the dispersion relation between po and the energy of 
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the particle. Indeed, for a free particle in the presence of a non-trivial chemical potential, 
the Green functions ()A.5jl and (jA.Tjl lead to 

S^ = i— — . +27r f{po,fi) 6 {p'^ -m'^), (3.6) 
p^ — + le 

where p^ = p'^p^ and, since the momentum integral of the first term vanishes, the charge 
density is obtained from the integration of the second term. Under these conditions, 
Eqs. ()3.3j) and ()3.4|1 are consistent with the definitions of the particle density n^^^ and 
current J(^) in (I2III). 

For fermions of opposite chirality (say left-handed ones) the previous expressions hold 
just changing ^ a^. On the other hand for several flavors, the mass eigenstates are an 
admixture of the weak eigenstates, with a mixing that depends on the value of the varying 
Higgs background. The expression of the diffusion equations and the sources, Eqs. ()2.2H) 
and ()2.22|1 . are related in a very simple way; namely the sources are trivially obtained 
by demanding a self-consistency relation in the absence of particle changing interactions. 
This property should hold, in first approximation, for any flavor structure of the theory. 
The relevant sources for the diffusion equations should be obtained by generalizing the 
trace in ()3.4j) to flavor space and eventually including a projection operator inside the 
trace, "P, when a given contribution is to be picked up. As stressed above, this formalism 
makes contact with the one proposed earlier by the authors of Ref. [H] and should lead to 
a more realistic estimate of the baryon asymmetry of the universe than the one presented 
in previous analyses of this subject. 

Although, as we have already outlined, we agree with the results of Ref. [H] in the 
formal definition of the sources, we shall now show that the structure of these sources 
differs from the one presented by the authors of Ref. ^T]. We shall try to clarify the 
origin of the discrepancy with those authors' results. 

For the chargino sector of the MSSM we define the right-handed and left-handed 
fermions as 

^-("^ = (l^) ' ^ ("^r ) ■ 

and expand the mass matrix to first order in derivatives around the point z as 

M(x) = M{z) + {x- zYM^{z) , (3.7) 
where we use the notation M^{z) = dM{z)/dz^ , and the mass matrix M{z) is defined in 

We now consider the first term in ()3.7p as part of the unperturbed Lagrangian, and 
the second term as the perturbation 

Cintix) = {x- zY [i,^^{x)M^{z)i,L{x) + V^l(x)Mt(z)V^^(x)} . (3.8) 

or, in the basis of mass eigenstates, 

Xr{,x) =U{z)ijji{x), xl{x) = V{z)ipL{x) 
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CU^) = {x - zY {x^^{x)U{z)M^{z)V\z)xl{x) + xi{x)V{z)Ml{z)U\z)xR{x)] (3.9) 



We denote by S^^ , S^'^, S^^ and S^^ the left-left, right-right, left-right and right-left 
Green functions of free fermions in the mass eigenstate basis with mass mi{z) as given 
in ()A.5|) . In this basis the right- and left-handed Green functions get modified to S^^ 
and S^^ by the presence of the interaction term ()3.9p and can be defined as solutions of 
Schwinger-Dyson equations. To first order in the insertion given by the interaction term 
fl3.9|l the solution to the Schwinger-Dyson equations is 

S^\x, y; z) =5««(x, z) + A5^«^(x, z) 

AS^^ix, y-z)= j d^w{w - zYS^^{x, w; z) U{z) M^{z)V\z) S^^{w, y; z) 

+ I d^w{w- zYS^^{x,w-z)V{z)Ml{z)U\z)S^^[w,y-z) (3.10) 



and a similar solution for . 

The CP-violating Green function l\S^^{x^y\ z) defined in ()3.1()j) can be written as, 

AS^^ = j [ ^p^e'P-'{a''GUM^M^U^G + cT,p''G''U{MMl + M^M^)U^G-h.c.} 

(3.11) 

where the notation = dG/dp^ has been used. 
Working out Eq. ()3.1ip it can be cast as, 

\qRR_\_ f ^^P Jp-r j ( G^l-^ll G1G2A12 

-4y(27r)4' Y''\-G,G,A',; -GlM[^ 



+ (a-p)(G'i^G2-G'iG'2^) ( ^([2^1 (3.12) 

The calculation of the coefficients J^x-, A12 and is straightforward. We will only give 
the relevant pieces that contribute to the CP violating current ()3.3|1 . 

Im [At,] 

Im [(M2U1 + fXU2)A'l2] 

Im [(M2M1 + fiU2) 

We can see from ()3.13j) that the combination U1U2 — U2U1 appears in the off-diagonal 
coefficients of AS^^. Had we neglected the latter we would have obtained only the 
orthogonal combination U1U2+U2U1 as claimed in Ref. However, contrary to the claim 
made by these authors, we do not find that the off-diagonal terms lead to a contribution of 
second-order to the currents relevant for electroweak baryogenesis. Since this is a crucial 
difference between our results and those obtained in Refs. [HI 112] , we will be more explicit 
about the origin of the discrepancy with those authors at the end of this section. 



Im M2/i ^ ^ 

[U1U2 + U2U'l) 

Im(M2/i) ^ ^ 
A -j^ [uiu'2 + M2M1 ) 

Im(M2/i) {U1U2 — U2U1) (3.13) 
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Returning to the derivation of the currents, the last step would be going from the 
Green function fl3.12j) to the CP- violating current as defined in (I3.4|l . To this end we 
must rotate to the weak eigenstate basis and project onto the Higgsino component. We 
must then define the Green function 

sI^jI = U^AS^^U (3.14) 

and the CP-violating background current ()3.4|) is now given by 

ilg'^ = limTrP22a^^(S(z + r/2,^-r/2) (3.15) 

r— >U 

where V22 = (o"o ~ '^3)/'^ with cij being the two by two Pauli matrices in flavor space and 
(Jq being the two by two identity matrix. The current jj^^ ^ can be given by the expression 

+ 2p^(Gi,G2-GiG2.) Im[(M2Mi+/iM2)i3i'^2] } (3-16) 

where the coefficients ^n, A12 and B12 of the Green function in ()3.13|1 are made explicit. 
The final expression can be written as, 

(M1M2 + ■U2M1) 



JRR 



Im(M2/i) f cf^p 



4A 



(27r)^ 



G\ - G\ - ^(Ci 



G2? 



+ 2p^ {Gi,G2-GiG2u) (niM^-n2<)} (3.17) 

The left-handed current j'^ij]^ can be obtained in the same way. The corresponding 
result can be read off from ()3.17|) after making the changes A ^ A and ui ^ U2. From 
them the vector and axial currents can be computed and the result obtained from ()3.17|) 
agrees with that presented in Eqs. (3.14) and (3.15) of Ref. |39j . 

To conclude this section we would like to make some final comments concerning the 
origin of our discrepancy with the results presented in Ref. |12j. In Ref. |12| it is correctly 
argued that, in the absence of particle-changing interactions and in the mass eigenstate 
basis at a given point z, terms proportional to WiWg — U2u'l only appear in the off-diagonal 
elements of the two-by-two matrix of Green functions. This is in agreement with the 
result obtained in Eq. ()3.12|1 . From here, the authors of Ref. jl2] conclude that the 
sources receive no contribution proportional to U1U2 — U2U1 at first order in the derivative 
expansion. This conclusion is not correct as we will now explain. 

First of all, let us emphasize that the particle changing interactions are essential in 
order to convert the original Higgsino density into a left-handed quark one. In the absence 
of these interactions with the plasma the net result for the baryon asymmetry will be 
much smaller than the one required for consistency with Big-Bang Nucleosynthesis. The 
Higgsino states can not be associated with the mass eigenstates at any point z along the 
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bubble wall, where the Higgs background is non-vanishing. In order to treat the particle 
changing interactions one has to transform to the weak interaction basis. 

Second, observe that the currents are completely determined by AS^^ (AS^^). The 
dominant, diagonal contribution 5''^^ (S^^) plays no role in the definition of the CP- 
violating densities at this order. The CP-violating densities receive contributions from 
both the diagonal and off-diagonal terms after rotating to the weak eigenstate basis as 
in fK114j) and both contributions appear at first order. Therefore, one should compute 
both the diagonal and off-diagonal contribution to the Green functions. The fact that 
the diagonal terms have no dependence on the combination — U1U2 does not imply 
that the sources share such an independence, since the off-diagonal terms in the mass 
eigenstate basis play an important role when interactions are included. Actually, as 
shown in Ref. [201, depending on the parameter space, the contributions proportional to 
U1U2 — U2u'^ may be the dominant ones in the generation of a relevant baryon asymmetry 
of the universe. 

Similar arguments can be used to show that the considerations of the authors of 
Ref. [Hj, leading to the absence of a dependence of the sources on M2M1 — can not 
be sustained. Indeed, Ref. jl2j tries to provide a formal proof to these considerations. 
Moreover, the authors of Ref. neglect the chargino mixing in the interaction terms, 
in spite of the fact that mixing is essential for a non-vanishing source. Therefore, their 
treatment relies on an approximation that is invalid in the interesting region where M2 
is close to /X. As we showed in Ref. ^^H] it is precisely in this region of parameters where 
the sources proportional to — U1U2 become most relevant in the baryon asymmetry 
generation (see discussion below). 



The detailed evaluation of the currents necessary to compute the sources for the diffusion 
equations presented in the last section was already provided in a previous article [39 . In 
Ref. the resulting baryon asymmetry was computed under the assumption that the 
connection between the sources for the diffusion equations and the CP-violating currents 
was that presented in Refs. [SSlEn], namely Si = rjn[^\ In section 2 we showed that, 
contrary to this assumption, the real sources, Eq. ()2.22j) . have a form that depends on 
derivatives of the CP- violating currents computed in the presence of the background fields. 
Therefore, the resulting value of the baryon asymmetry is different from the one computed 



The baryon asymmetry resulting from the solution of the diffusion equations depends 
on integral functions along the wall of the temperature dependent value of the sources 
times exponential factors which are slowly varying along the bubble wall fS^. Therefore, 
the new values can be simply obtained from the earlier ones by integration by parts, 
taking care of the proper boundary conditions. 

As discussed in the introduction, the total baryon asymmetry depends on two different 
contributions. The first one is proportional to 



4 Implications for the baryon asymmetry 



in Ref. 



eijHid^Hj = v\T)d^{l3). 
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This expression, Eq. M.lj) . is proportional to the variation of the angle 

13 = arctan [v2(T) / vi(T)] 

at the wall of the expanding bubble, which tends to zero for large values of m^- Further- 
more, independently of tua, for large values of tan/3, f3 varies only slightly from its value 
/3 ~ 7r/2 and therefore the baryon asymmetry tends to be suppressed. Finally, this con- 
tribution to the sources of the baryon asymmetry has a resonant behaviour for M2 = |/^|- 
Therefore, it becomes more relevant for moderate values of tua and tan /3 and for values 
of M2 ~ 

The second contribution to the baryon asymmetry depends on 

H^d^H2 + H^d^H^ = cos(2/?)a^/3 + vd^v sin(2/3) (4.2) 

Similarly to ()4.H1 . the first term in the expresion ()4.2|1 is suppressed for large values of 
niA and/or tan (3. Although the second term is also suppressed for large values of tan /3, it 
is unsuppressed for large values of m^. Therefore, we expect this contribution to become 
the dominant one for large values of m^, particulary in the non- resonant regions with 
values of M2 very different from 

Finally, let us remark that the final baryon asymmetry is obtained by solving the 
diffusion equation for the total baryon number. This equation may be easily derived by 
summing up the diffusion equations for the thirty six chiral quarks of the Standard Model 
(three generations, two flavors of quark per generation, two chiralities and three colors 
per quark). The sum of the quarks densities is just three times the total baryon number 
density (each quark carries a baryon number 1/3). The variation of the baryon number is 
governed by the sphaleron processes, which affect only the left-handed chiral quarks and 
leptons. 

The left-handed quark chemical potentials receive two contributions. The dominant 
one may be obtained by solving the diffusion equation for the different colors of quarks 
in the presence of gauge, Yukawa, mass and strong sphaleron interactions. Since no 
baryon number violating processes are included, the solutions to these diffusion equations 
lead to an equal number of baryons of a given chirality and antibaryons of the opposite 
chirality. These dominant densities may be considered as approximately constant during 
the characteristic long times in which the weak sphaleron processes take place. 

The left-handed densities receive also a subdominant contribution coming from the 
weak sphaleron interactions. This contribution is associated with a net-baryon number 
which, considering effective mixing between the different flavors and colors of quarks, is 
shared in approximately equal parts by all 36 of them. Finally, there is also a net lepton 
number created which due to charge conservation is equal to the net baryon number. 

Taking all these considerations into account, and considering the right-handed leptons 
out of equilibrium, one arrives to the equation 

/o ,,diff/ \ \ 

Dn"j,{z) - v^n'j,{z) = 9{-z)T^, i + AnB{z)\ (4.3) 

where F^^ = (a^ is the weak coupling constant), with k^s — 20 j4j being the 

weak sphaleron rate, while ix"^^ is equal to the sum over the three generations of the left- 
handed up and down quark chemical potentials associated with a given color of quarks. 



12 



as obtained by solving their diffusion equations in the absence of the slow sphaleron 
interactions. In the numerical simulations, we have considered heavy all squarks and 
sleptons except for the lightest stop, in which case the coefficient A = 24/7, which differs 
only slightly from the SM result A^^^ = 15/4 (see also Ref. p[Ti). 

In Fig. H we show the comparison of the baryon number to entropy ratio f] computed 
within the present model with the mean value of the one consistent with Big Bang Nu- 
cleosynthesis (BBN), TjBBN — (6 ± 3) X 10^^^ |ini. We have chosen soft supersymmetry 
breaking parameters in the stop sector such that they lead to a value of the Higgs bo- 
son mass consistent with the present experimental constraints, and the lightest stop light 
enough so that the phase transition becomes strongly first order, v{T)/T > 1 ^43] . 




100 200 300 400 



|i (GeV) 



Figure 1: Plot of r]/r]BBN a function of ^ for tan/? = 10, M2 = fi and the specified 
values ofrriA- 

Since the dominant component of the sources comes from the chargino sector |3()[ IHH] , 
the results of the baryon asymmetry depend only indirectly on the stop parameters, 
through the values of v{T)/T and I3{T). Different choices of the parameters lead to 
variations of the final result by a factor of order one. Since the method of computation 
has implicit uncertainties of similar order, so far v{T)/T > 1, the specific choice of the 
parameters does not affect the results for the baryon asymmetry in any significant way. 

We have chosen a maximal value of the phase of the fi parameter, siiKp^ = 1. There- 
fore, the inverse of the values shown in the figure can be interpreted as the value of sin (p^ 
necessary to obtain a prediction consistent with BBN. Observe that due to the uncertainty 
in the value of tibbn jlO], the value of tibbn niay be a factor 2 smaller than the above 
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quoted mean value and therefore the value of the phase may be a factor 2 smaller than 
the value obtained by the procedure described above. Hence, from Fig. [T] we see that 
for rriA consistent with the present experimental constraints, values of sin0^ > 0.05 are 
preferred. 



M2 = 200 GeV 




200 400 600 800 



[L (GeV) 



Figure 2: Plot of rj/rjBBN as a function of /i for M2 = 200 GeV, tan/? = 10, and the 
specified values of itia- 

In Fig. 121 we present the values of the ratio ti/tibbn as a function of /i, for a fixed 
value of M2 = 200 GeV. For small values of the dominant contribution comes from 
the resonant behavior at M2 = |yu|. However, as rriA increases the maximum of |?7| is 
obtained for values of n much larger than M2. The dominant, non-resonant component 
has opposite sign to the resonant one and for large values of fi compared to M2, it becomes 
only slightly dependent on rriA, in agreement with our discussion above. We see that even 
for large values of m^, values of the baryon asymmetry consistent with BBN predictions 
may be obtained for phases of order one. Observe that the baryon asymmetry tends to 
zero for large values of n, since the Higgsino component of the lightest chargino, present 
in the plasma, tends to zero in this limit. 

The above results are important in view of the constraints coming from electric dipole 
moments jl^I- It has been recently observed that, even in the presence of very heavy 
squarks and leptons, the electric dipole moment contributions may not be small, due to 
the existence of two loop diagrams involving chargino loops and the non-SM like Higgs 
bosons [46j. These contributions are enhanced with increasing tan/?, but become rapidly 
small for large values of rriA and small mixing in the chargino sector. The results of Fig.|21 
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show that a relevant baryon asymmetry may be obtained even in this particular regime 
of parameters. 

5 Conclusions 

In this article, we have derived the diffusion equations governing the result for the baryon 
density generated by the passage of the expanding wall at temperatures close to the 
electroweak phase transition temperature. Important for this result was the identification 
of the sources for the diffusion equations in terms of the currents generated by the non- 
trivial Higgs backgrounds, which vary along the bubble wall. We demonstrated that there 
are two different contributions to these sources, one proportional to e^^ Hid^Hj and the 
other proportional to Hid^H2 + (2 — 1). Contrary to previous claims, we showed that 
none of these contributions to the sources vanish. We clarified the origin of the discrepancy 
of our results with those of Refs. [HJ |121- 

The two contributions to the sources have a different dependence on the parameters 
of the minimal supersymmetric standard model. While the former contribution has a 
resonant behaviour for M2 = |/i| and goes to zero for large values of ttia, the latter has 
no resonant behaviour and remains no n- vanishing for large values of ttia- 

The final result for the baryon asymmetry shows that consistency with Big Bang Nu- 
cleosynthesis demands the relative phase between /i and the gaugino mass parameter M2 
(or Ml, if we consider the neutralino contribution) to be sizeable. Our calculations lead 
to a bound on the value of the CP-violating phase, sin0 > 0.05, but due to the natural 
uncertainties associated with our simplified treatment of the baryon asymmetry calcula- 
tion, and the experimental error in the determination of f]BBN, we cannot reliably rule 
out the possibility that somewhat smaller values of the phase 0^ may lead to consistency 
with the BBN predictions. 

On the other hand, for sizeable values of the CP-violating phase sinc/)^, cancellation of 
the potentially large electric dipole moment contributions is required. These contributions 
are dominated by one-loop effects with scalar particles and gauginos participating in the 
loop. These corrections may be efficiently suppressed for large values of the scalar masses 
of the first and second generation [121 • Under these conditions, non- vanishing corrections 
still appear at the two-loop level [IHI. The most relevant two-loop corrections are enhanced 
for large values of tan /?, but are suppressed for large values of the CP-odd mass parameter 
rriA- It is reassuring to observe that even for large values of the CP-odd mass parameter 
a baryon asymmetry consistent with BBN may be obtained. 

In summary within the framework of the minimal supersymmetric extension of the 
Standard Model, the realization of the electroweak baryogenesis scenario requires the 
presence of a light Higgs boson, with mass smaller than about 120 GeV, and a light stop, 
with mass smaller than the top quark mass |llj-ji28], Moreover it has been found 

that this baryogenesis mechanism does not guarantee a new CP violating signal at the 
5-factories Therefore, apart from the experimental constraints coming from electric 
dipole moments, a definitive test of the electroweak scenario withih the MSSM will come 
from Higgs and stop [4QJ searches at the Tevatron and at the LHC. 



15 



Acknowledgments 

We acknowledge discussions with T. Prokopec and S. Weinstock. M.Q. and M.S. would 
like to thank the Theory Division of the Fermi National Accelerator Laboratory, where 
part of this work was done, for hospitality. M.C. and C.W. also wish to thank the 
CERN Theory Division, where part of this work was done, for hospitahty. M.S. would 
like to thank the Theory Group of Institute de Estructura de la Materia, where this 
work was partly done, for hospitality. M.Q. would also like to thank the Physics and 
Astronomy Department of the Johns Hopkins University for the hospitality extended to 
him as Bearden Visiting Professor. The work of C.W. is supported in part by the US 
DOE, Div. of HEP, Contract W-31-109-ENG-38. Fermilab is operated by Universities 
Research Association Inc. under contract no. DE-AC02-76CH02000 with the DOE. The 
work of M.S. is supported by the US DOE Contract DE-A505-89ER40518. 



A The charging sector 

In this appendix we summarize some useful formulae which concern the chargino sector 
of the MSSM. The chargino mass matrix is given by 



y Ui{Z) He J 

where we have defined Ui{z) = gHi{z). The diagonalizing matrices are 



U 



A + A M2 Ml + /X* M2 



V'2A(A + A) V -{M2U1 + HeU2) A + A 



1 ( A + A M2U2 + llcUl\ . 

v/2A(ata) V -(^2«2+/i:«i) A+A ^-^ 

where field redefinitions have been made in order to make the Higgs vacuum expectation 
values, as well as the weak gaugino mass M2, real, 

/^^{Ml-\iJi,\^-u\ + ul)l2 
K^{Ml-\^x,\''-ul + u\)/2 

A = (A2 + |M2Xii + /x>2|')'^' , (A.3) 
and the mass eigenvalues are given by 

(A + A + »K^))-^^ + "iW"2(^K 



mi{z) 



v/(A + A)(A + A) 



^ . . ^ (A + A - Ul{z)) /Xe - U^{Z)U2{Z)M2 

V(A + A)(A + A) 
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In the mass eigenstate basis the free fermionic Green functions can be written in terms 
of the bosonic ones z) as 

S^'^ip; z) =a,p^G{p- z) S^\p- z) = (^"^f^ ^° G{p; z) 

S^^ip; z) = h*^'^ J,^ G{p; z) S^\p; z) =a,p^G{p; z) (A.5) 



and 



where Gi is the two-by-two matrix of a bosonic propagator corresponding to the mass 
|mj(z)| ()A.4j) given by 



= P^-fF{P^-Pr) 

Gf = [e{p') - h] {Pt - Pf) 
Gf = [0i-p') - If] {P.^ - Pr) 

Gf = -Pr - h {P^ - Pf) , (A.7) 

where fp = np^p'^l) is the Fermi-Dirac distribution function in equihbrium, which con- 
tains the dependence on the temperature T, 



' pl-p^-\mi{z)\'^±2iTi\p'^ 
and Fj is the particle width. 
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